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ABSTRACT 


The effect of gravity levels up to 400 times the local 
gravity value (a/g) was experimentally investigated in a 
Aero racing). cylindrical boiler. Heat flux levels up 
to 80,000 BTU/hr-ft* were reached. 

Observed improvement in heat transfer at low heat flux 
values with increasing gravity levels confirmed the results 
of previous investigations. At high heat flux values, a 
consistent trend was not observed. With increasing gravity 
levels, some data indicated a decrease in heat transfer 
while other data indicated an increase in heat transfer. 
Changes in heat transfer with variation of the liguid 


С о thickness was noted. 
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אד ה 


The study of the effect of increased gravity levels on 
nucleate pool boiling has been underway for approximately 
eleven years. This study has heen pursued because boilers 
and vapor-generation devices are being extensively revised 
because of new requirements for their use in spaceflight 
applications and in the widening use of nuclear reactors 
as energy sources. These boilers are required to operate 
at high heat flux levels and under variable gravity levels. 
The study of the effects of high gravity levels on the 
boiling phenomenon may also lead to a new understanding of 
the boiling heat transfer mechanism. 

Шис 15ם- 00 לוט 60רכ‎ investigation was to obtain im- 
5 data om nueteate boiling of water at high heat flux 
levels and at gravity levels up to 400 times the local 
Пру лао а/с). This was done in order to determine 
ІШ there was a Significant variation in heat transfer rates 
in this region with gravity level change and if there was 
amy GeverSal in the direction of the variation. 

Coste reena Tuthill I] studied the НЕО бї in- 
сасе Jraviryetevells of 1 to 45 a/g upon nucleate boiling 
heat transfer at heat flux levels between 100,000 and 
2007900 BTU/hr rA ПИ рала ЕЯ range of study, in- 
creased gravity levels required larger values of superheat 
for a constant heat flux. Heat transfer was decreased for 


imereases in acceleration in this region. 





Merte and Clark [2] studied the influence of increased 
acceleration, from 1 to 21 a/g, upon nucleate boiling heat 
transfer ın saturated water. Under pool boiling conditions, 
they detected an acceleration effect which was dependant 
upon the heat flux level. For values of heat flux less 
ан 20,000 BTU/hr а increased acceleration decreased 
the superheat required. At heat flux levels above 50,000 
BTU/hr a increased acceleration increased the superheat 
required. They postulated two different modes of heat 
meester tO explain this; nucleate boiling and natural 
convection heat transfer. At low levels of heat flux, 
natural convection heat transfer governed heat removal and 
increased gravity levels enhanced heat transfer rates. At 
high heat flux levels, nucleate boiling heat transfer was 
predominant and increased gravity levels decreased heat 
transfer rates. 

In the natural convection region Ref. [3] recommends 


Using formula 


үші 5 И 


Nu = C(GrPr : NUM ErPpr) < 2x10 


tor horizontal flat plates where C = .54, m= .25, Gr = 
Grashof number, Pr = Prandtl number, and Nu = Nusselt 
gBATX3 | | 
number.  Grashof number = 2 and =סען 0017 ב תו ה‎ 
V 


portional to the gravity level. This shows that as gravity 
level increases, heat transfer rates in the natural convec- 


tion region also increase. 





Adelberg [4] compared gravity forces with three other 
forces which also influence the heat transfer rate in 
nucleate boiling. These forces were bubble dynamic forces, 
Surface tension, and forced convection induced drag. He 
then reviewed existing literature for experimental data 
which he considered relevant and compared his analysis 
puru this data. 

Later, Adelberg [5] developed three equations for heat 
transfer in pool nucleate boiling. The range of applica- 
bility of each equation depended upon the gravity level. 
From experimental data available, he concluded that his 
equations for low and moderate a/g levels agree favorably 
with the limited data available. No data:was available to 
verify his equation for high a/g levels. 

Graham and Hendricks [6] investigated the influence of 
acceleration upon the mechanism of heat transfer in satu- 
Ес הט‎ ето from 1 to 10 a/g. They demonstrated 
that increasing acceleration tended to shift the boiling 
characteristic curve in the direction of improved heat 
transfer rates in proportion to the 1/4 power of the dif- 
ference in wall temperature and fluid bulk temperature. 

Gray, Marto, and Joslyn [7] investigated heat transfer 
behavior in a rotating boiler using acceleration up to 200 
ey aeene heat Gluxes up to 505,000 BTU/hr ft“, Ша СЕУ 
confirmed that et low heat fluxes, increased acceleration 
increased heat transfer rates; whereas at high heat fluxes, 


the opposite was true. This was in agreement with previous 





experimenters. The region where this crossover effect took 
Peace 00000760 at progressively higher heat flux levels for 
лек ассејека лопз. They concluded that unresolved ef- 
fects such as surface condition and aging affected the data 
at high heat fluxes as much as did gravity variations from 
005 200 8/0 . 

Marto and Gray [8] investigated high acceleration and 
heat flux effects on nucleate boiling. They achieved ac- 
celerations up to 400 a/g and heat fluxes to 818,000 BTU/hr 
ft^. Their results sbowed that for the inception of nucleate 
boiling from natural convection, heat flux and required wall 
superheat increased as acceleration increased. They also 
verified that at low heat fluxes, increased acceleration in- 
creased heat transfer rates. But at high 2. fluxes, heat 
transfer rates were practically independent of acceleration 
variations. In this region changes in heat transfer rates 
were as small in magnitude as the usual scatter in experi- 
mental data. 

Judd and Merte [9] studied the heat flux predictions of 
six different models proposed to explain the nucleate boiling 
phenomenon. They showed that none of the models adequately 
predict Hide СЕС ИОСаАШ О heat flux for other than saturated 
Pta 86ב 3075 3הבסס‎ standard gravity. All the models pre- 
dicted decreasing values of heat flux as acceleration in- 
creased. These observations suggested the omission of some 
mechanisms which are inoperative under standard conditions 
OL Gravity, but very sensitive to the level of subcooling 


and acceleration. 








Пете вето Han templirical procedure for correlating 


the results used in their study. Their equation for heat 
flux is 
. 1/3 ,98k3 1/3 Bar, ГІ. 
k(T __,-T,) пат e 
6 АІ АЕ УВ wall satv 


А Јо x 10 Pr — |C 
። (е те) 


ferric units must be used in this equation. 
unc רכה‎ 5 кетп e contribution of natural convection 
Enc tbe second term the contribution of nucleate boiling. 


They showed that as a/g increases, ከ ee decreases when 


both ШІ” ЕШ О/А ате constant. At low heat flux 


levels, the first term is important since natural convection 


dominates. As a/g increases, the heat flux may either in- 


сеток Честеазе depending upon the magnitude of ን.‏ )ה 


መ) At high heat flux levels, the 


compared to 
first term becomes insignificant and the second term domi- 
mates. In this region Q/A will decrease with increased a/g 
since ш un Tp) term decreases. This suggests that 
there may Бе а crossover region of heat flux where the ef- 


fect of increasing gravity level is reversed. 


l0 





PEO ከ! CONFIGURATION 


И оте с this High gravity nucleate boiling investi- 
Eon, the boiler section of an existing rotating thermo- 
syphon apparatus was used. 

Figure 1 is a schematic diagram of the thermosyphon 
assembly and associated equipment. שו רו‎ 2 75 a Cross 
sectional drawing of the thermosyphon. Figure 3 is a 
paurosraph of the eqgurpment setup in the laboratory. 

The eguipment was the same as that used by Woodard [10] 
with the following modifications. The entire thermosyphon 
was manufactured from pure copper stock with the exception 
of the two end plates of the evaporator. These were manu- 
Пәсіепгса from 304 Stainless steel. The grooves were designed 
to reduce heat transfer from the barrel of the evaporator 
пе Cnd plates. The difference in the thermal conduc- 
ורוו‎ O ОШОО GQitferent materials served the same pur- 
pose in this design. The condenser end was left open to the 
atmosphere by removing the .25 inch threaded plug in the end 


Section. 


A. EVAPORATOR MANUFACTURE 

The evaporator section was manufactured from pure copper 
stock and the end sections from 304 stainless steel. It had 
an inside diameter of 3.125 inches and an outside diameter 
Б ОСОБ :пецев. The length of the copper section was 3.125 
inches, and when assembled the overall length was 3.75 inches. 


The actual heated length was 3.5 inches. The heating coil 


Ug 





Was 11=gauge Chromel-A heater wire in a 0.1875 inch inconel 
ENMcath helrcally wound around the evaporator. 

The procedure for bonding the stainless steel end plates 
to the copper section and bonding the heater wire to the 
evaporator was an exacting process. Its purpose was to en- 
mera unitorm heat flux axially along the heater surface 
and to prevent electrical shorting between the heater wire 
and other components. The procedure was as follows: 

1. A thin sheet of gold copper foil was tack welded 

to the outside evaporator surface and each end 
section. 

2. The heater wire was wound over the gold-copper 

foil with gold-copper wire wound intermittently 


between the coils. 


3. The stainless steel end plates were then pressed 
places 
4. The entire assembly was brazed in a hydrogen 


atmospherically-controlled oven. 

5.  Sauereisen cement was air dried around the out- 

side of the heater wire. 

6. Asbestos insulation was wrapped around the heater 

with glass tape.‏ סי וסור וש 

Power was supplied to the heater coils by a 130 volt 
Dec. source Capable of delivering a maximum of 100 amperes. 
This power was transmitted by a brush assembly consisting 
of eight spring-loaded graphite brushes through two bronze 


rodeo the heater wire. The bronze collector 


ша 





калас ህር ር ራሙ supported by twelve 0.25 inch stud bolts. 
Phenolic spacers were used to align the rings and to in- 
late the bolts from the rings. 

The outside face of the evaporator was closed by two 
pyrex glass view ports. The inner port rested on a 
neoprene O-ring and was separated from the outer port by 
a compressed fiber gasket. A glass tube was fused to the 
memes port and passed through a hole in the outer port. 
This tube was used for evacuating and filling the system 
when used as a thermosyphon. 

In order to measure the temperature gradient in the 
evaporator wall, four holes, 0.0625 inches in diameter and 
2.375 inches deep, were drilled radially in the evaporator 
wall. Four sheathed thermocouples were inserted into these 
holes. The center of one hole was located 0.0625 inches 
from the inner wall and another hole's center was located 
the same distance from the outer wall. These holes were 
located 180 degrees apart. The remaining two holes were 
drilled so that their centers were 0.125 inches from the 
inner wall. These two holes ico located 180 degrees 
EEG. 

Due to the relative softness of copper, these holes were 
not drilled true, and it was thus necessary to measure the 
misalignment of each hole. This was done to determine the 
exact position of each thermocouple junction. The following 


procedure was used on each hole to determine its misalignment. 


ша 





₪ מנס 6565הה0‎ 2.5 inches in length was 
hardened at one end for approximately 0.5 inches to pre- 
clude any bending. This pin was inserted in the drilled 
hole with the hardened end outward. Protťtruding from the 
hole was a 0.125 inch section. A precision square 2.375 
inches high was aligned to the edge of the drilled hole 
ה‎ 305 the protruding pin end for alignment. When pro- 
perly aligned, the precision square was clamped in place 
and the metallic alignment pin removed. Then a Tungsten 
heliarc electrode approximately six inches in length and 
0.0625 inches in diameter was inserted into the hole. This 
electrode was chosen because it is manufactured to be pre- 
cision straight. The misalignment between the precision 
square and the tungsten electrode at the top of the block 
was measured using precision feeler gages. Since the 
electrode was assumed to be perfectly straight, the mis- 
alignment of the electrode at the top of the precision 
square was the same as at the bottom of the drilled hole. 

This procedure showed that the ends of the drilled 
holes were located as follows from the centerline of the 
оо SO 672, 1.683, апа 1.770 inches, respec- 
tively. The hole nearest the inner wall was only .016 
inches from the wall. A slight protrusion was visible. 
Due to the softness of the copper, the drilling may have 


бшсе е та Surface to bulge very slightly. 
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TNSTRUMENTATION‏ כ 
оте Measurements‏ ביר וו 

Temperature measurements were obtained from four, 
Ag auge, cCopper=constantan Е These were in- 
050 in 0.0625 inch diameter holes drilled 2.375 inches 
deep in the evaporator wall. They were sheathed in 0.0625 
memes OD stainless steel tubing. 

The manufacture of these thermocouples required a 
EDI procedure. First, the outer plastic insulation 
was removed from the wires by using a low flame to melt the 
insulation. The residue was then carefully scraped from 
the inner wires using a surgical blade. Extreme care was 
used so that the insulation of the inner wires was not 
ruptured. This outer insulation was removed from a one 
foot section of each thermocouple. 

Sie есе ессе! cubana, 0.063 inches OD, was cut 
into lengths of approximately five inches. The strippeo 
thermocouple wires were inserted into the tubing and the 
ends of the wires were soft soldered to make the junction. 
With this junction made, the wires were pulled back until 
only approximately 1/4 inch remained exposed. A coating 
of Conley-Weld Epo ከከ ን... placed over the exposed 
wires and junction. Then the wires were again pulled back 
צ1תס הנכהה‎ thetepoxy covered junction was exposed. The 
Sas Us cOmO protect the junction and to hold the 


wires securely in the tube. After the epoxy had set, fine 
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sandpaper was used to make a smooth finish and reduce the 
58 757 che tip tO that of the tubing. 
The leads from the four thermocouples in the evapor- 
ator wall were brought down the outside of the condenser 
and through the rear condenser flange to a junction board 
at the end of the shaft. From the junction board, the leads 
went inside the main shaft to a slip ring assembly which 
was attached to the main shaft by a coupling. 
The output from the slip rings was fed to a Hewlett- 
Packard 20106 Data Acquisition System. This system provided 
ес ок visual output of the thermocouple readings. 
የጠ. | .፡!፡፡!።.] Speed Measurement 
The rotational speed of the evaporator was measured 
Шипа а magnetic pickup and a Systron Donner Counter 1013 
Series. An electronic strobe light was used as a check 
ое ооо Water Temperature Measurement 
הנפמ הבה ה ו‎ cooling water temperatures to 
the condenser section were measured using a Hewlett-Packard 
2801A Quartz Thermometer. The uncertainty in these readings 
мае +.1°6. 
4. Ns ets Power Measurement 
The power into the evaporator section was measured 
using DC voltage and amperage meters. 
5. Cooling Water Flow Rate Measurement 
The cooling water flow rate to the condenser was 
Measured using a calibrated flow rotameter. The calibration 
procedure for this Е Ча described in бет. ПТ. Тһе 


uncertainty in these readings was +.5%. 
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6. Liquid Annulus Thickness Measurement 
The liquid annulus thickness was measured by a paper 
rule scribed in 1/16 inch marks attached to the viewing glass 
at the evaporator end. The uncertainty in these readings 


. פסתסתג 5.032 = 


PROCEDURE‏ שא[ תתל 

The evaporator was filled with approximately 200 ml of 
distilled water using a syringe. Filling was accomplished 
either from the condenser end or the evaporator end thru the 
glass tube. When using the latter method, the glass tube 
was sealed after filling. Cooling water flow was then es- 
kablıshed to a 30% reading on the rotameter. The drive 
motor was then energized and the desired RPM established. 
When a liquid annulus had formed, DC power was applied to 
the evaporator heater. With the desired RPM and power level 
established, data was recorded when thermal equilibrium was 
established. Waiting time between data recording ranged 
from 5 to 15 minutes after a change in variables. The fol- 
lowing data were recorded: RPM, a/g, power supplied to the 
пеасег, tour thermocouple readings, coolant mass flow rate, 
temperature difference across the condenser, atmospheric 
pressure, and liquid annulus thickness. 

Securing the machine was accomplished by first shutting 
off power to the evaporator heater and then shutting off 


the drive motor. Lastly, the cooling water was secured. 
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Pie eR MENT AT, RESULTS 


пегпой OF heat flux calculation was to use‏ | הר רו 
the four thermocouples installed in the evaporator wall to‏ 
determine the radial temperature distribution across the‏ 
evaporator wall. With this temperature distribution deter-‏ 


mined, the inside wall temperature, T could be computed 


wall’ 
along with the heat flux, 0/2, 85 described in Appendix B. 

Power level was limited to 6 kilowatts. Above this 
level, steaming occurred from the condenser end and con- 
siderable vapor was lost to the atmosphere. The initial 
data run revealed that with the highest power level avail- 
able, the temperature variation across the evaporator was 
too small to be accurately determined by the thermocouples. 
This was due to the high thermal conductivity of the copper 
evaporator, the thin wall, and the limit on the power level. 
The errors in the location of the thermocouple wells de- 
СООЛ у Section Tk also contributed to the uncertainty 
in the thermocouple readings. Therefore this method of 
calculation was abandoned. 

The alternate method used to calculate the evaporator 
wall temperature and the heat flux was to make a heat 
balance on the condenser section. The heat supplied to 
the evaporator equals the voltage times amperage supplied 
to the heater. The heat to the condenser section equals 
the heat supplied to evaporator minus convection and radia- 
tive heat losses to the atmosphere. The heat removed from 


the condenser section equals the heat supplied from the 
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evaporator section plus the heat generated by friction in 
B wears and bearings plus the heat generated by viscous 
dissipation in the coolant. These last two sources of heat 
generation can be combined into the heat generated due to 
rotation. An overall heat balance states that the heat 
supplied to the evaporator equals the heat removed from 

the condenser minus the heat generated due to rotation. 

The measurement of the heat due to rotation was deter- 
mined at three different RPM's. Cooling water flow was 
established at a constant rate and electric power to the 
evaporator heater was secured. When steady. state 
conditions were established, the temperature difference of 
the cooling water and the mass flow rate were recorded. 

The heat removed from the condenser was calculated using 
the mass flow rate and temperature difference across the 
condenser as described in Appendix B. 

Evaporator wall temperature, wall saturation temperature, 
and heat flux vere then calculated by the second method de- 
scribed in Appendix B. The specific data taken were enumer- 
ated under "Operating ыы All reduced data appear 
in Appendix E. 

Data runs 1 thru 5 were taken at constant RPM with vary- 


PMigevahucs OQumneat tlux, Баба taking began at 25 G's up to 


400 G's. The liquid annulus thickness decreased from .4 
בכרה‎ ОС ШІ |І ЕО 219 inches at run 5. Data runs 6 thru ll 
were also taken at constant RPM and varying neat flux. In 


ріне ሪክ) 6 started at 400 G's and decreased to 25 G's 
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during run 11. The liquid annulus thickness decreased from 
«ከ ሮ ከይር during run 6 to .38 inches during run 11. Data 
шише 12 Enea 16 were taken at constant heat flux with vary- 
ing values of RPM. Liquid annulus thickness decreased from 
.38 inches to .25 inches during this series of runs. Data 
Ai was taken at a constant gravity level, 25 G's, and a 
Semstent heat flux of 80,000 BTU/hr 4 This run was taken 
to determine the effect of liquid annulus thickness of the 
heat transfer coefficient. Liquid annulus thickness started 
at .5 inches and decreased to .125 inches. Data run 18 was 
taken at 400 G's with varying values of heat flux. Liquid 
annulus thickness started at .63 inches and decreased to 

.56 inches. This run was taken to determine if there was 
any effect upon the data due to increasing the waiting time 
after each variable change to at least 15 minutes. This was 
well beyond the time that steady state appeared to have been 
established from the indicator values. 

ПАСЕ ШЗ פמסנ%הלהבבמס‎ of the state of the liquid in the 
evaporator were made during each data run. A stroboscope 
instrument was used to light the evaporator interior. At 
low gravity Levels agitated nucleate boiling took place in 
the liquid annulus at all heat flux levels. The degree of 
בשו‎ a d aa heat flux increased. At high gravity 
levels the nucleate boiling process increased very slowly as 
heat flux levels were increased. The degree of turbulence 


was much less at high gravity levels. Any effect of liquid 
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. Р 
nnulus thickness variation upon nucleate boiling was not 


шоу арракеп from visual observations 
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PA ISS TON OF RESULTS 


Figure 4 shows the results of data runs 2 thru 5. RPM 
was held constant and heat flux was varied. Gravity levels 
started at 50 G's and progressed to 400 G's. Liquid annulus 
thickness was approximately .2 inches during this series. 
This figure shows that at low heat fluxes heat transter is 
improved as gravity levels increase. There is no crossover 
region discernible, although the data appears to merge to- 
gether at high heat fluxes. 

Figure 5 shows the results of data runs 6 thru 9. RPM 
was again held constant and heat flux was varied. Liquid 
annulus thickness was approximately .5 inches. Data а 
started at 400 G's and decreased to 50 G's. Again higher 
gravity levels show improved heat transfer rates at low 
values of heat flux. In this figure there is a definite 
035 the 100, 200, and 400 G curves. The data at 
50 G's appears to be offset at superheats which are too 
mich. 

Бајине ломе the results of.data runs li thru 16. 
Liguid annulus thickness was approximately .38 inches during 
these runs. Heat flux was held constant and RPM was varied. 
There is a definite improvement in heat transfer with in- 
creasing gravity level, but no crossover region is evident. 

Пао пета ስዕ effect of liquid annulus change on 
the heat transfer Coefficient. As the thickness of the 


annulus decreases, the heat transfer coefficient increases. 
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This agrees with the results of the investigation of MacKenzie 
ғала the results of Ref. 13. 

Figure 8 illustrates the effect of extending the waiting 
РФ ПЕ аг сасп Gata point from 5 - 7 minutes up to 15 minutes. 
This figure is a comparison between data runs 6 and 18. It 
appears that the longer the waiting time the more reliable 
the data. Even though the parameter indicators show that 
steady state has been established after approximately 5 
Minutes, waiting until at least 15 minutes results in a 
larger temperature difference and a higher heat flux. Thus 
it takes much longer than is readily apparent to reach 
Бају state conditions. 

Figure 9 is a comparison between data collected at 25 G's 
during these tests and data presented in Ref. 2 at 20 G's 
and in Ref. 8 at 25 G's. There is good correlation between 
0 0 5005 02 Кес. 2, but the data of Ref. 8 has higher 
values of superheat. 

0 33 פסַהַהּע"5‎ 56 08068 collected at 50 G's with 
that presented in Ref. 8 at the same gravity level. The 
ወ. кош Век. 8 15 offset in the direction of higher 
superheat. A eebe explaination for this is different 
surface conditions and ace. berweenstheir test 
surface and the copper surface used here. 

Figure 11 is a comparison between the 400 G data during 
this test and the 400 G data presented in Ref. 8. There 


is good correlation between both investigations. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


CONCLUSIONS 


Е. 


Puras Mmueleate Pool boiling at low heat flux 
levels, there is a definite improvement in heat 
transfer with increasing gravity levels. 

At higher heat flux levels some data shows a 
crossover trend while other data appear to just 
merge. No definite conclusion can be drawn from 
the available information. 

Increasing the waiting time after any variable 
change beyond the apparent time of steady state 
enhances the reliability of the data. 
Variations in the liguid annulus thickness 


Е ን ፕ ውን ን) ኤስ ሮመ תג‎ the heat transfer coefficient. 


RECOMMENDATIONS 


| 


Increase the waiting time after any variable 
change to at least 15 minutes to increase the 

ми шара су об сре data. 

Operate with the machine closed to the atmo- 
sphere so that steaming is prevented and heat 
flux levels can be raised above those obtained 
in chas investigation. This will necessitate 

a device for measuring the temperature or pres- 


спе evaporator vapor space.‏ וה כ 
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Provide fins on the outside of the condenser to 
increase the heat transfer rate to the cooling 
Weser, 

Increase 6 evaporator wall thickness so that 
larger temperature differences Can be realized. 
Install a permanent photographic installation 
so that the boiling phenomenon can be visually 


recorded. 
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APPENDIX A 


CALIBRATION OF THERMOCOUPLES 


The copper-constant thermocouples were calibrated using 
the Hewlett-Packard Data Acquisition System, 20T0C, for 
Шо: readout. The following fixed point calibrations 
were made: 

Ace POInt of water (32.0 °F) 

гот па point of water (212.0 °F at 14.696 psia) 

ЕИ Тесто ропе оғ їп (449.4 °F) 

Melting point of lead (621.1 °F) 

Calibration of the ice point of water was accomplished 
by using a Thermos bottle filled with crushed ice and dis- 
tilled water. Each thermocouple was inserted in this ice 
bath and three separate sets of data were recorded. The 
average value of these three data sets was used. 

For boiling point of water, a large test tube was filled 
with distilled water, and a loose fitting cork with a hole 
in the center was inserted. Each thermocouple was inserted 
through the hole in the cork to approximately 1/2 inch 
above the water level when violent boiling had occurred. 
Four separate sets of data were taken and the average value 
used. An uncertainty for each thermocouple was also cal- 
culated using the highest and lowest reading. The largest 
value was £.009 millivolts. The standard temperature was 
corrected for atmospheric pressure variation using the 


= ei onstormula Given in Ref. [7]. 
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Maita cion tac the melting point of tin was accomplished 
by heating a pure tin sample in a large test tube. A small 
diameter glass tube with one end closed was inserted into 
the sample. Then the thermocouple to be calibrated was 
placed in the glass tube. A high range glass thermometer 
was also inserted in the sample. The temperature of the 
sample was then raised above the melting point and the heat 
withdrawn. The molten sample was allowed to cool until the 
freezing temperature was reached. Then for a discernible 
amount of time the temperature stabilized. At this point 
data was taken. This procedure was repeated six times 
morder to insure accurate results. The average of the 
six data points was used and an uncertainty value was then 
Calculated for each thermocouple. The largest uncertainty 
value was +.004 millivolts. 

The same physical setup was used to calibrate the thermo- 
couples at the melting point of lead as was used for tin. 

0 005 1688 was cooled in the molten state to its melting 
point, it did not stablize for any apescctable time. . There- 
ШО И рат сен was obtained using the minimum time 
interval possible. This data was then plotted on graph 
paper and the area of minimum slope was used as the calibra- 
tion point. The average reading in this area of minimum 
slope was determined and uncertainty value again Calculated. 


The largest uncertainty value was +.1 millivolts. 
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Reference 14 recommends a third order equation of the 


OEM 


i, = E ae bt? ns 


for interpolating thermocouple readings between 100 ?C and 
approximately 300 *C. Using this equation and the calibra- 
tion temperatures, the constants were determined using the 
IBM 360 computer for each thermocouple. Then using this 
result, a table of thermocouple readings versus degrees 
Centigrade and Fahrenheit vas calculated for each Teme 
couple using the computer. 

The results of the interpolating equation were plotted 
every 50 °C inorder to compare these values with the ex- 
perimental data. At the three experimental points the un- 
certainty values were also plotted. The calculated results 


compared well with the experimental data. 


28 





КЕРИМОЯЛХ- В 


DATA REDUCTION 


DC ALCULATION OF EVAPORATOR WALL SATURATION PRESSURE 


Evaporator 
a oo. ፦”።፦ (Беше ле 


liquid annulus 
h 


К LTDA TIT P Wall 


evaporator inner radius 
n distance from centerline to liquid annulus 
ћ - liquid annulus thickness 
Therefore 
room re > h 


both sides of the equation‏ 00028 ה 


ДЕ 2 = г. - ar sh + m 


Dnezrassumptions in this calculation were that solid 
ped rotation takes place and that the density of the two 
phase mixture could be approximated by the density of water 
ו‎ е Ot rotation of a fluid moving as a 
ቢት ያ፡፡ ያ פגכב פה‎ is Called forced-vortex motion [15]. 
The pressure field varies with radial and axial position. 
Since the machine was maintained in the horizontal position 
шеге мав по ахла! variation in pressure, and since rota- 
tion was at high RPM, earth's gravity was neglected. The 


radial component of the pressure becomes: 
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= = (1) 
ox gc 
where Py = Сеше רה‎ 010 at 212 °F, Lbm/£t* 
w = angular velocity of evaporator, rad/sec 
и racial position, ft 


ПОШ Шо ме стоп tactor, 32.17 ООН ርን ኔም 


P = pressure, БОГ ТС” 


INECgrating equation (1) from га to Y: yields: 








р. м 
o (r.° = x“) 
1. а 
ос 
Eubstituting for 5. yields: 
2 
up 
ДР = ———h (2r. - h) 
дас = 
dr 
2 
от“ 
Nc i im т 2 к к 
а 2gc 
2 - 2 
where P = pressure at wall, Lbf/ft 
satw 
Es = pressure in vapor space, ШОО ОО 


B. CALCULATION OF EVAPORATOR WALL SATURATION TEMPERATURE 
The subroutine for the method of least squares on the 
IBM 360 computer was used to find an appropriate polynomial 
וק‎ 205 tne water Saturation pressure-temperature curve. 
)וו‎ ОЕ 14.0 to 50.0 psia were used and a fourth 
order polynomial gave excellent correlation between actual 


data and a number of selected test points. 
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The polynomial used was as follows: 


ПИ 50070: = со. о2б42Р | - .168663P |, ^ 
satw satw satw 
+ .0024795P__, 9 - .00001488665 ^ 
ሽ satw ` satw 
where 
0 = saturation temperature at evaporator wall, °F. 


ACA LCULATION OF TEMPERATURE PROFILE 

The heat flow through the evaporator wall during steady 
Beate Conditions can be reduced to a one dimensional problem 
05 1S nO axial or circumferential heat conduction. 
Under these assumptions the steady state conduction equation 


in Cylindrical coordinates is 


.א 


1 d e. 
= A 2 
where r = radial coordinate 
Te temperature 
ወ ስሞነ ን ን በጋሪ Conductivity of copper. 


integrating equation (1) once gives 


ar 2 E 
בע (2) א‎ = constant = В, 


The thermal conductivity is assumed to vary linearly 


with temperature. Therefore, 


Bene = го Е м(т-то) (2) 


> 





Data on the conductivity of copper was obtained from 
КОШ 16 . In the temperature range of interest, O0 ?F to 
КУ Ене conductivity of copper is a linear function of 
temperature. This data was plotted and ko and m, the 


slope, were determined graphically. 


ko = 231.77 BTU/hr-ft-°F 
m = -.02222 оре 
ШӘ - 0 ФЕ 


ENESEItutring equatron (2) into (l) and then integrating 


again gives: 


eon fame 2-To)JT = By + B, ה ה‎ 


The subroutine for the method of least squares on the 
IBM 360 computer was used to find the constants В, апа Во: 
Using these constants and value of the inner wall radius, 


Т was extrapolated. 


O A CULATION OF HEAT FLUX, Q/A 


Prom Fourier's Law: 


| ат 
КА ПТ к ay! 


0 


Or 


27114 Ва 


ርጋ 
| 


The inside area of the evaporator is 


А = 2T"r.L 
I. 
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Then 


d Q l PI 
A 2тг. Г. ie ! hr £t2 


ШИРСАТСЦОГАТТОКМ OF HEAT FLUX, Q/A, (SECOND METHOD) 


ШЕ ሰባን тег rate out of the condenser is 


out = п ל‎ и ШЕ = 
where m = coolant mass flow rate EN) 
Ct, = coolant specific heat oe) 
T4 = coolant inlet temperature (°F) 
T, поене пне temperature (°F) 
S: computed was then corrected for the heat loss due 


Eomevaporator rotation. 


NC corrected 


pac orehedrted surface 


Then Q/A for the evaporator - 


Dan corrected 


zone t, 
1 


B. CALCULATION OF EVAPORATOR WALL TEMPERATURE, = 


- 


(SECOND METHOD) 


For steady state conditions, the heat flow equation is 


У T ) 
One - wall [17] 
1n _ © 
Igo 
L 
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where Ша = average of the two center thermocouple reading, 


КЕ 


1 


15 tadial distance tO thermocouples, 


С 


ו 


Then dividing by the inside area of the evaporator gives: 


зт ) 
ug 
SE С - wa 
ההש‎ = 
ү TN 
il 
пошту па for и yields: 


Le 
= (1n = (Q/A) 
i 


Tvali © te T К 


mine value of К was assumed constant at 226. 
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ПЕВТ Г 





С‏ או 


UNCERTAINTY ANALYSIS 


In order to determine the reliability of the data col- 
lected and to determine the most probable error, an uncer- 
tainty analysis was performed utilizing the Kline and 
McClintock technique [19]. The analysis was performed on 


ner data taken during run 1 at 1 KW. 
PRO NCERTAINTY IN THE CALCULATION OF HEAT FLUX 


п m cp AT 


„ [w2 w 2 ка? W,2 
Wq’ _ m ep AT A 


W = uncertainty in subscripted quantity 
2 





10,045 BTU/hr-ft 
eon 1,005 BTU/Lbm °F 
a. ft" 


m = 305 


И = 1,5% 

War E 

W = 1 
Cp 

Wa = 1 


Substituting into equation (1) yields: 


= (197 = 4.875 
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IN THE CALCUALTION OF EVAPORATOR WALL‏ וו 


TEMPERATURE 





= „130 1 ft 

1 
226. / BTU/hr-ft-?Fr 
ul = +.001 
r = .1398 ft 

С 
Wrc = t2 

= + 
W Br 
Go 10,045 BTU/hr-ft“ 
= 1 

Te 
Wo, = 220487 

q 


substituting into equation (2) yields: 
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UNCERTAINTY IN THE CALCULATION OF WALL SATURATION 


PRESSURE 


= P 
satw atm 29 





2 


גלש 2 





W = / (М יי‎ 


(2r,-h)W_) 
Psatw Patm 29, 4 p 


ПО Ln —— 











2 
- A 2 р V? ж-ж 7 2 
А (2r; h)W 2) + 2422 В (25, п) ИЗ.) 


ПИ TTT 





2 2 
how 2 рм 2 
ы (--- М.) ( бл 111.) (>) 
( ge ri 26 1 h 
ВИ. 1302 ft 
1 


М = 78.5 rad/sec 


Wg, = +.001 

h - .0338 BTU/hr-ft^-ep 
Woop, = 23.0 

Нұ, = 1 

ЕЕ 221 ft/sec 

W, = 5.05 

W, = +.032 

L = 60 NECS 

3 = 2 א 


ስህ unto equation (3) yields:‏ רו 


Wp 
БЕНЕН ዕጩወያታ ኤሽ. ሽል -25% 
satw 


D 








ШЕБЕТАТУЦТУҮ ТЫ ТПЕ CALCULATION OF WALL SATURATION 


TEMPERATURE 
T О ИСЛ 160592012 ғ ን ከባ: ረ 
satw satw satw 
4 .0024795 P 5 = 0000148866 P 4 
у са ту : багы 


(4) 
Substituting into equation (4) yields: 


WTsatw Ж 8 


O‏ ו 


=A 015 = O 


The largest uncertainty appears in the calculation of 


ШОО са: flux. This is due to the error in the mass flow 


rate measurement and to the error in the coolant flow tem- 


perature difference measurement. 
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